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We recently reviewed the experimental determination of the volume fraction, 4>, of hard-sphere 
colloids, and concluded that the absolute value of <j> was unlikely to be known to better than ±3- 
6%. Here, in a second part to that review, we survey effects due to softness in the interparticle 
potential, which necessitates the use of an effective volume fraction. We review current experi- 
mental systems, and conclude that the one that most closely approximates hard spheres remains 
polymethylmethacrylate spheres sterically stabilised by polyhydroxystearic acid 'hairs'. For these 
particles their effective hard sphere diameter is around f-fO% larger than the core diameter, de- 
pending on the particle size. We argue that for larger colloids suitable for confocal microscopy, the 
effect of electrostatic charge cannot be neglected, so that mapping to hard spheres must be treated 
with caution. 



I. INTRODUCTION 



A collection of hard spheres is one of the simplest ex- 
amples of an interacting system. Hard-sphere packings 
have been important since the dawn of civilisation [l], Q , 
while using hard spheres to model the liquid state dates 
back at least to Kirkwood jH in the 1940s. The sim- 
plicity of hard spheres lends itself to analytical theory 
[3-01 and computer simulation, pl-fioj While the absence 
of inter-particle attraction precludes a liquid phase, con- 
centrated hard sphere fluids capture many properties of 
the liquid state. [I|,[ll],[l2| Due to its analytic tractability, 
and that it is described by just one parameter, volume 
fraction <j), the importance of the hard sphere system as 
a basic model in understanding condensed matter cannot 
be overstated. 

Early hard sphere experiments include Bernal's use 
of ball bearings to model liquid structure. 01 However, 
ball bearings have negligible thermal motion. Suspen- 
sions of mesoscopic colloids exhibit Brownian motion and 
are thermodynamically equivalent to atoms and small 
molecules [l4j . In 1986, Puscy and van Mcgen (l5j 
showed that a suspension of sterically-stabiliscd poly- 
methylmethacrylate (PMMA) particles showed hard- 
sphere-like equilibrium phase behaviour. Their work led 
to many experimental studies of the statistical physics of 
hard spheres using colloids as models. Since Pusey and 
van Megen's work, the equation of state of hard-sphere 
colloids has been determined fl(| . crystal nucleation has 
been observed [U HH > an d the glass transition has been 
studied [H. 

The body of experimental research just reviewed relied 
on light scattering as the structural and dynamical probe. 
The advent of single particle tracking in real space with 



confocal microscopy j20t |21| opened a new dimension in 
experiments on hard-sphere-like systems, yielding an un- 
precedented level of detailed information [lH- Confocal 
microscopy of hard-sphere-like suspensions is thus ideal 
for studying generic processes where local events are im- 
portant, such as crystal nucleation (23|, melting Q and 
dynamical heterogeneity (2f| [26[ . 

In principle, the thermodynamics of a system of hard 
spheres is controlled solely by the state variable <j>. We 
have recently reviewed the experimental determination of 
volume fraction[23] and concluded that, although relative 
values of <j) may be known with high precision, absolute 
values can only be determined to within 3-6% accuracy. 
This matters, especially when dealing with dynamical 
properties (e.g., phase transition kinetics and the glass 
transition), since these can be very strong functions of <j>. 

But the accurate determination of <f> is only part of the 
experimental challenge. The other part of the challenge 
was hinted at by the title of Pusey and van Megen's 1986 
paper, "Phase behaviour of concentrated suspensions of 
nearly hard colloidal spheres" [l5j , where we have added 
the italics to emphasise the point in question, namely, 
that true hard spheres do not exist in reality. In Puscy 
and van Megen's case, the lack of hardness is almost cer- 
tainly due to the small but finite compressibility of the 
PHSA stabilising 'hairs'. The same is generically true of 
other sterically-stabilised particle systems. 

In very-nearly-hard, sterically-stabilised suspensions, a 
new sample has to be made for every state point (p. Apart 
from being cumbersome, this also restricts the accuracy 
with which the sole thermodynamic control parameter 
can be 'tuned'. Thus, there has been a drive to use par- 
ticles such as 'microgcls', whose diameter is temperature 
dependent. Since temperature, T, can be tuned far more 
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FIG. 1: Schematic representation of various models for hard- 
sphere colloids, (a) Sterically-stabilized particle, with surface 
'hairs' (not to scale), where the average hair thickness 5 and 
the core-shell diameter a ca = o c + 28 are needed for a full 
characterisation, (b) Microgel particle, which is a heavily 
cross-linked polymer, (c) Charged colloid, where the electrical 
double layer (shaded) gives rise to an effective diameter a c g. 



accurately than (j>, this allows the scanning of a single 
sample very finely through (f> space by varying T. But 
the price paid for the 'tunability' of particle diameter is 
that some softness is built in by design. 

To understand a final reason why real model hard- 
sphere colloids may be somewhat soft, consider why in 
practice, not only <p but the particle diameter u mat- 
ters. Experimentally, colloids as synthesised are usually 
not matched in density to the dispersing medium. Signifi- 
cant sedimentation (or, less usually, creaming) during the 
timescale of an experiment therefore presents a problem. 
In light scattering studies, this is circumvented to a large 
extent by the use of small particles, say, a < 400nm, so 
that for PMMA particles dispersed in cis-decalin, a single 
particle sediments < 1mm per day. The middle part of 
a bulk sample would therefore be little affected by sed- 
imentation over a day. But if one wants to eliminate 
the effect all together (e.g. for long-time measurements), 
or when larger particles (say, a > lOOOnm) are needed 
for imaging purposes, then solvent mixtures for density 
matching are needed, which often introduces significant 
charging. The practical need for larger particles therefore 
inevitably brings softness. 

In this second part of our review, we survey critically 
these three sources of softness in experimental systems of 
colloids that have been used to model hard spheres. Fig. 
[TJ In the spirit of the first part of our review [231, we 
seek to provide a means by which the hardness may be 
assessed, and consider the consequences this may have 
on the behaviour of the system. Likewise, we argue for 
more clarity on the part of experimentalists, concerning 
the softness of the systems they use. We also suggest a 
number of criteria by which experimentalists (and theo- 
rists using experimental data) may judge whether a cer- 
tain system of colloids may be considered 'hard enough' 
for answering particular physics questions. 

Below, we first discuss mapping experimental systems 
onto hard spheres via measuring the inter-particle po- 
tential in section |TTJ We then treat softness of non- 
electrostatic origins in section [TTTI and of electrostatic ori- 
gins in section IIV[ in general, both effects are present 



simultaneously. Finally we give a worked example in sec- 
tion [V] illustrating many of the points raised throughout 
this review, before concluding in section [VTl 

II. MAPPING TO EFFECTIVE HARD 
SPHERES VIA THE INTERACTION 
POTENTIAL 

Given that real colloids inevitably display a degree of 
softness, it is important to be able to map their behaviour 
to that of hard spheres for the purpose of comparison 
with theories and simulations of perfect hard spheres. By 
mapping, we mean finding an effective hard-sphere diam- 
eter cr ff so that one may map from the experimentally- 
controllable particle number density p to an effective 
hard-sphere volume fraction 4> c s using 4> c s = TTpa^ s /6. 
There are two conceptually distinct ways of determin- 
ing <7 e ff. First, one could map via some known hard- 
sphere property, such as the volume fraction at freez- 
ing (4>f = 0.494), or the viscosity of the suspension as 
a function of 4>. This class of methods have been re- 
viewed in detail in the first part of our review. (27| Here 
we concentrate on a second class of methods: determin- 
ing (7 e fT from the inter-particle potential, u(r). To do so, 
of course, requires means of measuring it(r), which is the 
main topic of this section. 

First, however, we briefly review how a knowledge of 
u(r) can be used to determine an effective hard-sphere 
diameter, cr e ff. Perhaps the simplest approach is to set 
an effective hard sphere diameter akT such that the inter- 
particle repulsive energy at this centre-to-centre separa- 
tion between two particles is equal to the thermal energy, 
i.e. 

0u(r = a kT ) = 1, (1) 

where j3 = 1/ksT. A more sophisticated approach, 
which distinguishes between different potentials with the 
same akT, is to use the Barker-Henderson effective hard 
sphere diameter [28| 

oo 

obh = Jdr[l- exp (-/3u(r))} . (2) 
o 

Other, yet more sophisticated mappings, exist, such as 
that due to Andersen, Weeks, and Chandler [2{j]; this 
approximation is known to work well for mapping the 
static properties of liquids and more recently structural 
relaxation time near the glass transition [3(3 • All of these 
approaches rely on knowing u(r). We now review meth- 
ods for gaining this knowledge. 

A. Direct measurement 

A host of sophisticated techniques are now available 
for direct measurement of colloidal forces; integration of 
the measured force-distance relationship then gives u(r). 
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In the case of PMMA particles sterically stablised by 
poly-12-hydroxyl steric acid (PHSA) 'hairs', the inter- 
action may be inferred via the direct measurement of 
the interaction potential between mica surfaces coated 
by PHSA using the surface-force apparatus. The results 
were well described by an inverse power law, suggesting 
a reasonably (but not absolutely) hard interaction 3l|; 
see section UTTl for more details. 

Other methods for measuring colloidal interactions di- 
rectly include total internal reflection microscopy, which 
measures the force between a colloid and a glass wall,[32j 
and atomic force microscopy with a colloid on the tip 
of the cantilever. (33j The interaction between two non- 
index-matched colloids confined to a line can be measured 
by optical tweezers. [lij 

One attraction of such direct methods is that no a 
priori assumption need be made about u(r). However, 
their use requires care. Thus, e.g., in the case of optical 
tweezers, relatively small and subtle experimental errors 
can lead to the wrong sign of the interaction between 
charged colloids (35| . 



B. Extraction from correlation functions 

For a colloid at volume fraction 0, the inter-particle 
potential, u(r), uniquely determines the system's pair 
correlation function, g(r), or, equivalcntly, the structure 
factor 5(q), which is essentially the Fourier transform 
of g(r). Determining g(r) or S(q) from a given u(r) is, 
of course, one of the classical problems of liquid state 
theory. [3g In principle, it is also possible to reverse this 
procedure, and infer u(r) from measured correlation func- 
tions. The inversion of S(q) to obtain u(r) has a long 
history in liquid state physics, and has also been used 
for colloids. [13, A fundamental difficulty with this 
approach is that, like all inverse problems, this one is ill- 
conditioned. Essentially, many different forms of inter- 
action can give rise to the same S(q). For our purposes, 
it is instructive to bear in mind that the S(q) for the 
inert gases near their respective triple temperatures can 
be well fitted to the hard-sphere S(q) at an appropriate 
4>, even though the inter-particular potential under the 
same conditions is well approximated by a Lennard- Jones 
form: u(r) = Ar~ 12 + Br~ 6 , which includes an attractive 
part, and a repulsive part that is very far from 'hard'. 

Inverting the real-space pair correlation function, g(r), 
can also give the interaction potential, [H, [4(j and is sub- 
ject to the same a mbig uities, especially when many-body 
effects are present. 4l| But if g(r) can be measured in the 
limit of vanishing </>, then the inversion to give u(r) is 
unique: [15 



lim g(r) = exp [ 
4>— >0 



-Mr)] 



(3) 



This approach has become possible with the advent of 
real-space techniques, which allows the determination of 
g(r) by direct counting. [2l], [4(| E3 In the limit repre- 
sented by Eq.O monodispersc hard spheres give a perfect 
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FIG. 2: Attractions in "hard" spheres for confocal mi- 
croscopy. Experimental data are from [4^]. The dashed black 
line is from Percus Yevick theory for hard spheres at <f> = 0.071 
[3^ |. The solid line is computer simulation data with parti- 
cle tracking errors and polydispersity added, for a square well 
attraction of depth fcsT and range 0.09u. 



step function form for g(r), but polydispersity and par- 
ticle tracking errors would blur the sharpness of the edge 
at r = cr, as would a small degree of softness. In Fig. 
[5J we show the g(r) measured in this limit for a puta- 
tive hard-sphere suspension. [43l| It is immediately clear 
that these particles cannot, in fact, be hard spheres. The 
peak in the dilute-limit g(r) is due to a short-range inter- 
particle attraction. [86j Computer simulations can be used 
to fit the measured dilute-limit g(r) using a square-well 
attraction for u(r); importantly, however, a distribution 
of particle sizes as well as random tracking errors (both 
modelled by Gaussians) are essential to obtain a good 
fit. Interestingly, it has been arguedQ that Eq. ((3|) re- 
mains a remarkably accurate approximation at finite but 
modest d>. 



C. Charged particles 

The application of Eq. [3] makes no a priori assump- 
tion about the form of u(r), and can therefore be applied 
to any system. For charged systems, if the inter-particle 
potential is modelled using, e.g. DLVO theory (for which 
see section IPV|) . then one could access u(r) via measuring 
the parameters of this theory, viz., an effective charge 
on the particles, Z, and the ionic strength of the sol- 
vent. The definition and therefore determination of Z is 
highly non-trivial, [46| and the deduction of ionic strength 
from conductivity measurements in non-aqueous solvents 



depends on a number of assumptions 
reasonable results can be obtained 140 




Nevertheless, 
45 



D. Accuracy 

All of these methods have limitations. Those based on 
particle tracking methods (optical tweezers and extrac- 
tion from real space) are limited by the accuracy with 
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FIG. 3: Estimation of effective colloid-colloid interactions 
in sterically-stabilised PMMA particles: (a) oh = 200 nm, 
(b) oh = 2000 nm. In both parts, light pink lines denote 
u s (r), the interaction due to the steric stabilisation. The 
dashed blue line in each case represents unscreened weak 
electrostatic interactions in low dielectric constant solvents 
(cis-decalin/TCE) calculated for effective charges Z = 2 and 
Z — 16 for (a) and (b) respectively and a Debye length 
of k _1 = 5000 nm. The solid blue line represents typi- 
cal screened electrostatic interactions: (a) in water (charge 
Z = 1700 and Debye length = 4 nm), and (b) a density 
matching mixture of cis decalin and CXB (charge Z — 500 
and Debye length n^ 1 — 100 nm). 



which the particle coordinates can be tracked, typically 
30-100 nm, or up to 0.05cr for a 2 /im particle. More- 
over, there is evidence that tracking errors on two closely 
approaching colloids may tend more to the point of con- 
tact, rather than being uniformly distributed. [3ol |43| Fur- 
thermore, averaging over many particles in an inevitably 
polydisperse sample effectively adds a further error to 
the measurements. [43[ Comparison with simulations can 
be used to compensate for some of these errors, though 
this avenue has been little explored to date. Finally, at 
higher concentrations. (40j and in reciprocal space, one 
needs to make prior assumptions about the form of the 
inter-particle potential that one wants to determine in 
the first place. [36T - l38j Accurate measurement of u(r) is 
therefore far from straightforward. 



III. INTRINSIC SOFTNESS 

We now turn to review softness in model hard-sphere 
systems of non-electrostatic origin. There are two generic 
sources of such 'intrinsic' softness: steric stabilisation, 
and the use of particles with temperature-tuneable size. 



A. Steric stablisation 

A sterically-stabilised particle is shown schematically 
in Fig. [TJa). We have previously [27] reviewed ways 
to arrive at an effective hard sphere diameter for these 
particles via mapping to various hard-sphere properties. 
As already mentioned in section III Al such mapping is 
also possible via direct measurement of colloidal forces. 



This has been done for PHSA-stabiliscd PMMA parti- 
cles; the measurements are quite well describe d by an 
inverse power law potential with energy scale e c :[3l| 



(t 



(4) 



The relative range of u s depends on the particle size, for 
example n = 170 for particles with a hydrodynamic di- 
ameter of an = 200 nm particles, and increases with par- 
ticle diameter. Likewise the strength of the interaction 
also depends on the particle size, with u s (an) = 146 fcgT 
for a h 



are 



200 nm. The results of Bryant et al. [31 
replottcd in Fig. [3] These quantify what is intuitively ob- 
vious, namely, that for a fixed length of stabilising 'hairs', 
larger particles are relatively harder. 



B. Microgels 

One main motivation to use soft particles to model 
hard spheres is that the diameter, and therefore concen- 
tration, of these colloids may be temperature tuneable. 
The most popular systems here are microgcl particles. 
As cross-linked polymer coils, their swelling is related to 
solvent quality, with the precise degree of swelling con- 
trolled by the cross-linking density. There are two va- 
rieties: 'neutral' microgels, in which a low amount of 
charges is screened by added salt or is so small as to be 
negligible, and ionic microgels, which carry (much more) 
electrostatic charge, for which see section ITVl 

Few analytic expressions for microgel-microgel interac- 
tions exist. Expressions for polymer-covered fiat surfaces 
have been suggested as a first approximation [49j . The 
steric interactions between neutral microgels have been 
likened to crosslinked polymers, so that a form for the 
interaction can be obtained provided the density profile 
of the particle is known. {50l| Additional interactions in 
ionic microgels do have analytic forms (5l| , which can be 
added to the steric repulsion, for which an inverse power 
form is often assumed. 

Irrespective of the absence of well-attested analyti- 
cal forms for the inter-particle repulsion, swellable col- 
loidal microgels are definitely not hard, and the degree 
of swelling depends sensitively on the experimental con- 
ditions, including possibly the concentration of microgel 
particles. If the latter is a significant effect, then any c/> e ff 
would become state-dependent. 

We now review a popular system, dispersions 
of poly r/V-isopropylacrylamid) (PNiPAM) microgel 
particles, [52| for which water is a poor solvent at T > 
33°C, so that at and above this temperature, the diam- 
eter of the particles in an aqueous environment dramat- 
ically shrinks. Changing the amount of the cross-linker 
iV,7V'-methylenbisacrylamide tunes how much shrinkage 
occurs. [1U Sufficiently monodisperse samples crystallise 
at high volume fractions, giving crystals whose struc- 
ture has been variously reported as face-centred cubic [54| 
or as consisting of the more or less random stacking of 
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hexagonally-packed layers. [Hj Richtering and co-workers 
have examined the physical properties of PNiPAM mi- 
crogel suspensions with a variety of probes and discussed 
their findings in terms of possible mapping onto hard 
spheres. We mention three aspects. 

First, neutron scattering shows that an individual par- 
ticle has a constant-density core, surrounded by a corona 
in which the density gradually decreases to zero over a 
distance that is approximately twice that of the core. [56j 
From this finding alone, we expect the particles to be 
significantly soft. 

Next, the structure factors of PNiPAM suspensions 
at progressively higher particle concentrations have been 
measured, and compared to that of hard spheres. (57| Us- 
ing two different methods of data analysis, it was con- 
cluded that the S(q) of these suspensions could be de- 
scribed adequately within a hard-sphere framework by 
assigning <^> e ff (or, equivalently, a c s) to the particles at 
concentrations <p e s < 0.35. Above this concentration, 
increasingly large deviations from hard-sphcre-like be- 
haviour were observed. Interestingly, in a later paper. [58| 
the same group points out that the mapping to hard- 
spheres at (/> e ff was obtained by making one of two as- 
sumptions: either that a hard-sphere structure factor 
(from Percus-Yevick theory) in fact fitted the data, or 
that the form factor of a single particle determined 
in the low-concentration limit did not change when the 
microgel concentration was increased, which is not self- 
evidently true. 

Finally, Richtering and coworkers investigated the 
fluid-crystal coexistence gap.f53T.l5ll In one study of PNi- 
PAM particles dispersed in water. [53j a 4>e« was deter- 
mined by requiring agreement with the hard-sphere ex- 
pression for suspension viscosity at low concentrations, 
r l/ 7 lo = 1 + 2.50 o ff + 5.9<^ ff (where 770 is the solvent vis- 
cosity). This procedure is problematic because of the 
aforementioned state-dependence of a c g (and therefore 
of (j> e ff). Nevertheless, using this mapping, the fluid-solid 
coexistence gap was found to be <// ff = 0.59 < <j) < 0.61 = 
4>^g, i-e. it occurs at significantly higher concentrations 
than that in hard spheres (0.494 = cj)^ < (p^ s < <p m = 
0.545), and is substantially narrower (<f>^ — (// ff — 0.02, 
</>gg — (/)/j S = 0.051). Comparison with simulations f59| 
of the freezing of particles interacting via a power-law re- 
pulsion u(r) oc r~ n gives n sa 13 for PNiPAM particles 
with 240nm < <th 300nm. This is considerably softer 
than the inter-particle potential found [3l[ for sterically- 
stabiliscd particle of comparable size (er c = 200nm), 
which can be characterised by a power law with expo- 
nent n = 170. 

A different mapping was used in a second study of 
fluid-crystal coexistence, 58] in which the the microgel 
particles are now dispersed in dimethylformamide (DMF, 
a good solvent chosen for refractive index matching). The 
freezing concentration of the microgel particles expressed 
in mass fraction, /i = m m i cl . ogc i/(m micrO g i + TO so i vcnt ), 
was converted to an effective hard-sphere volume frac- 



tion by a multiplicative factor, <p c R = Sfi, where the 
'swelling ratio' S was chosen to yield a freezing volume 
fraction of </){ ff = 0.494. Interestingly, this procedure 
gave 4>^ff ss 0.55 for the point at which 100% crystallisa- 
tion should occur, consistent with hard-sphere behaviour. 
This agrees with the findings by another group of a re- 
cent imaging study using lar ger PNiPAM particles dis- 
persed in aqueous medium, [55j but contrasting strik- 
ingly with the previous finding by the same group of 
a significantly narrower coexistence gap.[53[ The same 
study found that the collective diffusion of these micro- 
gel particles dispersed in DMF and their hydrodynamic 
interactions could not be well described by any mapping 
to hard spheres. 

These studies illustrate some of the difficulties associ- 
ated with mapping microgcls to hard spheres. An ad- 
ditional issue is that of polydispersity. Both softness in 
the inter-particle potential and polydispersity (60j af- 
fect the width of the coexistence gap, so that the effect of 
these two quite distinct physical factors may be difficult 
to disentangle. Fortunately, microgels can be synthesised 
with polydispersities as low as ~ 1% before swelling, so 
that perhaps the polydispersity effect can be neglected 
in the first approximation (cf. the very small effect 1% 
polydispersity on the miscibility gap of hard spheres (60|). 

A variant of the 'canonical' microgel has been synthe- 
sised and characterised by Ballauff and co-workers [6ll l62j 
consisting of a hard polystyrene core onto which is grafted 
a network of cross-linked PNiPAM, so that the swollen 
shell has approximately the same dimensions as the core 
radius (~ 50 nm). To determine </) c ff, a core volume frac- 
tion <f) c was first measured by conversion from the particle 
mass fraction using the density of polystyrene. The hy- 
drodynamic diameter of the particles, cth was then de- 
termined using dynamic light scattering; this was later 
confirmed to be very close to the diameter of the outer 
corona visible in cryo-transmission electron microscopy 
(cryo-TEM) images, [13] which also gave the core diame- 
ter cr c , and the polydispersity of the core-shell diameter 
distribution (w 9%). Finally, using (j> c {f = 4> c (o~h /o~c) 3 1 
the fluid-crystal coexistence gap at 21°C was found at 
0.483 ± 0.007 < cff < 0.546 ± 0.007, which is, within 
experimental uncertainties, very close to the hard-sphere 
interval of 0.494 to 0.545. 

It is interesting to analyse these findings further. As 
the authors themselves [62l| have pointed out, if <p^ s is 
rescaled to exactly 0.494, then melting occurs at 0.556. 
This gives a coexistence gap wider than in perfect hard 
spheres. If this does not reflect experimental errors, then 
the situation is somewhat unusual - the most common 
'culprits', polydispersity and softness in the repulsive po- 
tential, both narrow the coexistence gap [Fig. U (a) and 
(b)]. However, short-range attraction in the potential 
has the opposite effect of widening the coexistence gap if 
the polydispersity [H| is low enough (63l - l65| . On the other 
hand, 9% polydispersc hard spheres should be at or be- 
yond the experimental limit of crystallisation, [66j proba- 
bly due to the onset of multiple solid phase coexistence 
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in the phase diagram. [67j which requires long-range par- 
ticle motion for fractionation. That crystallisation was 
still observed in these PNiPAM samples to give a coexis- 
tence gap wider than that of hard-spheres underlines the 
the lack of complete understanding of microgel physics. 

More recently, Ballauff and co-workers studied in de- 
tail the rheology of a 17% polydisperse suspension of 
their core-shell particles at and near the glass transition 
(found to occur at <fi e ff = 0.640), and compared their data 
with mode-coupling theory (MCT) calculations for hard 
spheres. [6^1 At this polydispersity crystallisation was in- 
hibited, but this introduced an extra level of complexity 
into calculating <p c g. The authors relied on the fact that 
the high-frequency viscosity had been found to be rela- 
tively insensitive to polvdispersitv. (69j and mapped the 
high polydisperse system onto a less polydisperse system 
in which (f> e g had already been calibrated according to 
the procedure explained above. A large measurement of 
agreement with MCT predictions has been found. We 
note in this connection that the comparison with MCT 
mostly relies on a relative measure of the distance to the 



glass transition 



<j)g)/<pgi so that the work 



is perhaps less vulnerable to systematic or statistical un- 
certainties in arriving at </> e ff. 



C. Star Polymers 

In the spectrum of hard to soft colloids, microgels oc- 
cupy a middle niche. Towards the soft end of this spec- 
trum, we encounter star polymers; when the number of 
'arms' (or 'functionality, /) is high enough, these can 
be considered soft colloids, with decreasing softness as / 
increases. Thus, for example, stars with / > 34 are pre- 
dicted to crystallise [70|. But the structures of the crys- 
tal phases differ from the crystals nucleated from hard 
spheres (random stacking of hexagonal layers). More- 
over, the effective potential between two particles at 
close approach is logarithmic. (7l| which is significantly 
softer than the power-law forms encountered in this re- 
view so far. Thus, although in a mixture of (large) hard 
colloids and (smaller) star polymers, 32-arm stars were 
found to behave almost hard-sphere like vis-a-vis the 
hard colloids, [?H it is unlikely that pure star polymer 
suspensions can be used as models of hard spheres. 



The treatment of electrostatic interactions on the 
mean-field, or linearised Poisson-Boltzmann (PB), level 
is largely adequate for our purposes. This is especially 
true for non-aqueous systems, [JJ [zl| in which the ener- 
getic penalty of ionisation is high, so that ion densities 
are low, and multivalent ions can be safely neglected. 

The linearised PB theory is incorporated into the 
Derjaguin-Landau-Verwey-Overbeek (DLVO) theory 74 
to describe the interaction between charged colloids. The 
original DLVO potential consists of van der Waals (vdW) 
and electrostatic components. We are primarily inter- 
ested in situations in which sterically-stabilised particles 
become charged, so we will assume that the steric repul- 
sion is adequate to render the vdW component negligible. 
Instead, we will consider an inter-particle potential con- 
sisting of a steric repulsion, u s (r), and an electrostatic 
interaction, which in linearised PB theory has an Yukawa 
form, uy(r): 



u(r) = u s (r) + u Y {r), 

exp[-«(r - cr c )] 



u Y {r) = e Y - 



r/o c 



Here, the contact potential given by 



(1 + no c /2) 2 a c 



(5) 
(6) 

(7) 



where, Z is the colloid charge, and the inverse Dcbyc 
screening length is given by k=^/47tAb/9,: „, where pi on 
is the number density of monovalent ions. The Bjerrum 
length 



A; 



/3e 2 /(47re er), 



(8) 



is the distance at which the interaction energy between 
two electronic charges is fc_eT, where e is the electronic 
charge, eo the permittivity of free space, and e r the dielec- 
tric constant. While this form of the electrostatic interac- 
tion is only valid in the range that linearised PB theory 
holds (weak electrostatic interactions), higher charging 
can also be treated with a Yukawa interaction by using 
a renormalised charge that is smaller than the physical 
charge on the particles. [75j 



A. Strong charging in aqueous solvents 



IV. ELECTROSTATICS 

It is now generally accepted that immersion of a col- 
loid in a liquid medium always gives rise to some degree 
of charging of the particles. This is a potential source of 
softness that should always be considered in experiments. 
If charged groups are present on the surface of the parti- 
cles, entropy inevitably favours at least a certain degree 
of dissociation; even when charged surface groups are ab- 
sent, adsorption of charged species from the medium will 
give rise to charged particles. 



Water has an unusually large zero-frequency dielectric 
constant, e r w 80 at around room temperature, giving a 
Bjerrum length of 0.7nm. Since Xb is comparable to the 
size of small ions, many ionic salts are readily soluble in 
water. This in turn means that the Debyc length is eas- 
ily controlled down to a few nanometres. For a = 400nm 
and Z = 1700 in water with a = 4.0nm, Eqs. [51 E] 
and Eqs. Q]and[2]give out = 1.052<r c and obh = 1.062<7 C . 
Thus, screened charged colloids in water very often ap- 
proximate hard spheres reasonably if they are either ster- 
ically stabilised or the short k _1 is still large enough to 
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FIG. 4: (a) The theoretical phase diagram of hard spheres at 
different polydispersities, a. F = fluid, S = (crystalline) solid; 
thus FSS denotes fluid-solid-solid coexistence. Replotted from 
Wilding and Sollich l60|j. (b) Phase diagram of hard-core 
Yukawa particles, from [73] in the absolute volume fraction - 
Debye length (1/kct) plane. Here the contact potential ey = 
8fcsT. In the case of zero Debye length, the hard sphere limit 
is recovered. Replotted from Hynninen and Dijkstra [77l ]. 



render vdW attractions irrelevant. Nevertheless, some 
aspects of the behaviour of such colloids may still be 
non-hard-sphere-like. Thus, Piazza et aZ.[lj| found that, 
although the equation of state (EOS) of their screened 
charged colloids approximated the hard sphere EOS at 
moderate to high 0, at low 0, the barometric distribu- 
tion in the earth's gravitational field, <f>(z) = 0(O)e -z / z °, 
was not observed, with the measured height distribution 
showing a more slowly decaying 'tail'. This is related to a 
decoupling of ions and colloids, leading to a macroscopic 
electric field 1761. 



Weakly-charged small particles in non-aqueous 
media 



In non-aqueous solvents, charging is profoundly al- 
tered. The reduced dielectric constant, e < 10, leads 
to Xb ~ 10 — 40 nm. With ionic sizes in the range of 
1 nm, one expects strong coupling between oppositely- 
charged ions and and therefore little dissociation of sur- 
face groups. It was therefore a long-held assumption that 
electrostatics could be safely neglected. However, the 
work of the van Blaaderen group (zH and others have 
shown that some electrostatic charging always occurs. 

It seems that the degree of charging in many 
systems 0, H [zl Izl H3| can be described by the rule of 
thumb ZX bI o~ ~ 6. Thus, the system used by Pusey and 
van Mcgen.[15j sterically-stabiliscd PMMA (a 700nm) 
in a mixture of cis-decalin and carbon disulphidc (e r = 
2.64, Xb ~ 20nm), can be expected to be charged to 
some extent, as later work on a similar suspension seems 
to confirm. [801] 

Nevertheless, when the particles are small enough, such 
charging can often be ignored. Figure [3{a) (dashed red 
line) shows uy (r) for a charge of Z = 2 on the surface of 
<jh = 200 nm particles (corresponding to ZXb/o « 6) in 
a dispersion medium with a Debye length of k~ 1 = 5/xm. 
The measured steric repulsion for sterically-stabiliscd 



PMMA particles of this size[3l| is also shown. It is 
clear that for all relevant length scales in this situation, 
uy(r) <C ksT. These particles can plausibly be consid- 
ered hard spheres. 



C. Particles for confocal imaging 

The most popular hard-sphere model system to date 
for confocal microscopy is sterically-stabilised PMMA, 
because of the possibility of using a solvent mixture 
for simultaneous refractive index and density matching. 
For accurate determination of coordinates, particles with 
°~h ^ 1/im are required. The density-matching solvents 
used are typically halogenated hydrocarbons. The use of 
these solvents has a number of undesirable side effects. 
They tend to swell the particles much more aggressively 
than non-halogenated hydrocarbons, and they can dam- 
age the fluorescent dye molecules included in the particles 
for laser confocal microscopy. They also lead to signifi- 
cant leves of electrostatic charging. 

Four commonly used halogenated solvents are cy- 
cloheptyl bromide (CHB), cyclohexyl bromide (CXB), 
tetrachloroethylene (TCE), and carbon tetrachloride. [89j 
In these low dielectric constant solvents (e.g., e r = 7.9 for 
CXB), colloids suitable for confocal microscopy (ctjj ~ 2 
/im) acquire a charge of Z ~ 100—500. The Debye length 
in a density-matching CXB-cis decalin mixture can run 
to microns as the ionic strength (due predominantly to 
solvent self-dissociation) can be as low as 10~ 10 M, [78| 
which is much lower than the ionic concentration in pure 
water (10 -7 M). 

The charge on the colloid, although much lower than 
what can be expected on similar sized particles in water, 
is now almost unscreened, which can lead to very long- 
ranged and strong interactions (ey >100 fcsT). These 
interactions can and do vary from sample to sample, as 
the ionic strength in CXB and CHB varies from batch 
to batch, and as a function of time. 0, EE IHl__Coi- 
loidal crystallization has been found in some cases [23j at 
(j> w 0.4 (as compared to hard spheres at c/)/j S = 0.494), 
but at least one experiment saw crystallization at vol- 
ume fractions as low as <p ~ 0.01;[78| often these crys- 
tals were body-centered-cubic (bec) in contrast to hard- 
sphere crystals which arc random-hcxagonal-close-packcd 
(rhep) or fee. Furthermore, for some batches of CXB, 
the colloid charge can change with particle concentra- 
tion, leading to strongly (^ -dep endent interactions, and 
even to re-entrant melting. [47| The use of PMMA par- 
ticles in halogenated solvents for confocal microscopy is 
therefore problematic. 

The situation can be improved somewhat by the use 
of salts to screen the charges. The problem with this ap- 
proach is that salts soluble in these solvent mixtures such 
as tetrabutyl ammonium bromide (TBAB)[78| are solu- 
ble only to around 260 nM.j47l.l48j This results in a Debye 
length of k" 1 w 100 nm. Although this is substantially 
less the diameter of imageable colloids (<jh ^ 1/^m), it 
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is not negligible and a noticeable degree of softness will 
likely result, Fig. Eljb). On the other hand, since the ma- 
jority of ions now come from the salt, the ionic strength 
and therefore the colloid-colloid effective interactions will 
likely be reasonably independent of <j>. 

Another possibility is to use lower dielectric constant 
solvents such as TCE and CCI4 (e r = 2.5 and 2.24 respec- 
tively) . Lower e r increases Xb , Eq. [5J The rule of thumb 
for estimating the degree of charging, viz., ZAb/<7 c ~ 6, 
therefore predicts a lower Z. However, both of these sol- 
vents are strongly absorbed by PMMA, and can lead to a 
volume swelling of > 40%. [45[ Unless the swelling is very 
closely monitored and characterised, it becomes a serious 
source of potentially large systematic errors, [27j because 
(f) oc u . Solvent absorption also increases the density and 
refractive index of the particles. Thus, one of the initial 
attractions of using such halogenated solvents is lost - 
without swelling, adding one of these solvents can den- 
sity match nearly exactly but also (fortuitously) nearly 
match the refractive index. With significant absorption, 
more TCE (say) than is needed for index matching has 
to be added to achieve density matching. Unless a third 
solvent is used to re-achieve index matching (which itself 
may need to further swelling), a turbid sample results. 

Such turbidity not only degrades image quality, but 
can also give rise to significant vdW attraction. For 
example, Fig. (2[ which shows the measured g(r) of a 
(f) = 0.071 suspension of sterically-stabilised PMMA par- 
ticles in a density-matching mixture of cis-decalin and 
TCE. [4^] The pronounced peak at touching immediately 
alerts us to the presence of inter-particle attraction, as 
simulations confirm. 



V. A CAUTIONARY WORKED EXAMPLE 

Since a well-attested analytical form for the interaction 
between charged colloids is available, Eqs. ([S])-©, we 
close by presenting an 'exact' comparison between the 
simulated phase behaviour of such particles, Fig. [H and 
various ways of mapping their behaviour to hard spheres. 
This comparison illustrates how careful one must be in 
drawing conclusions from such mapping, even in the limit 
when inevitable experimental uncertainties (e.g. due to 
polydispersity or slightly soft steric-stabilising 'hairs') are 
negligible. 

Consider charged hard spheres, so that u s (r) in Eq. ([5]) 
is the perfect hard-sphere repulsion. We model scenar- 
ios that may reasonably represent sterically-stabilised 
PMMA used in confocal imaging, and take a = 2 /xm, 
and a Dcbye screening length of = 100 nm (so 
that k<7 = 20). The latter is a round figure chosen to 
correspond roughly to CXB with the maximum possi- 
ble amount of dissolved TBAB (260 nM). Consider two 
colloid charges Z = 500 and Z = 100. The former is 
consistent with the rule of thumb, ZA^/cr w 6, while 
the latter is 5 times lower than predicted by this empiri- 
cal relation. These charges have been reported in differ- 
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FIG. 5: Mapping the phase behaviour of two hypothetical 
monodisperse charged hard sphere colloids (with parameters 
based on real systems - see text for details) to pure hard 
spheres. The hypothetical particles and solvent have the fol- 
lowing properties: a = 2/Ltm, k" 1 = lOOnm. The particles 
have two different charges (a) Z = 500, (b) Z — 100. In 
each case, differently shaded and delimited regions denote the 
fluid-solid coexistence gap of pure hard spheres (grey, 'HS'), 
and from: simulations [73] (red, 'sim'), mapping using Eq. J5J 
(lilac, '<Tfcr')i an d mapping using Eq. ^ (blue, 'ctsh')- In 
(b), the result of mapping using Eq. <(Sj) does not change the 
coexistence region from that predicted by 'sim', and is not 
shown separately. 



ent studies [40J, |47J, |81|, [82j for nominally identical PMMA 
particles used for confocal imaging dispersed in a density- 
matching mixture of cis-decalin and CXB. Equation ([B]) 
gives ey ~ lOfc^T and ps 0.5fceT for these two charges 
respectively. 

Taking these parameters, the simulation results [7^ re- 
plotted in Fig. 0Jb) can be used to determine freezing 
and melting for our hypothetical systems. These are de- 
limited by red lines in Fig. [5] for (a) Z = 500 and (b) 
Z = 100. Note first that these values of not-very-large 
surface charge, variously reported in the literature for 
nominally very similar PMMA colloids, in fact pertain 
to rather large differences in the freezing/melting transi- 
tions, both in terms of absolute values of the transitions 
points and in terms of width of the coexistence region. 
Thus, some kind of 'mapping' is clearly necessary if we 
are to use either system to model hard spheres meaning- 
fully. 

We proceed to calculate cr c s for the particles repre- 
sented by these two parameter sets using either Eq. (fTJ) 
or Eq. ([21), and Eqs. flSJ-flSJl. If mapping is to be helpful, 
we should expect that once we have transformed <p to 
</> e ff, freezing and melting of the two system should occur 
at approximately the hard-sphere values, viz., 0.494 and 
0.545. 

Fig. [S{a) shows that for the case Z = 500, Eq. (fTJ) pre- 
dicts freezing at <f>* kT = 0.484, just 0.010 from the 'correct' 
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value. [9C§ However, the width of the coexistence gap is 
reduced, which is the result of the softness of the screened 
electrostatic (Yukawa) repulsion. Since the mapping in- 
volves scaling </> by constant, </> c ff = (<7 c ff/er) 3 0, it pre- 
serves the relative coexistence gap, so that the melt- 
ing concentration, (j) kT is significantly underestimated. 
Conversely, the Barker-Henderson treatment gives a rel- 
atively accurate estimate of (j)™ T , but rather significantly 
overestimates </>{ T . In the case of Z = 100, since 
€y < kgT, OkT = so that Eq. ([1} maps perfectly 
onto hard spheres. In this case, the Barker-Henderson 
approach does the worse job, even when compared to the 
raw (unmapped) coexistence gap given by simulations of 
the bare Yukawa interaction. We note that since the co- 
existence gap varies with the interaction, mapping to the 
true hard sphere volume fraction at freezing inevitably 
gives an erroneous melting volume fraction unless the 
colloids are absolutely hard. 

Figure [5] therefore shows that the different methods of 
mapping to hard spheres do not give the same result. 
In practice, of course, the inter-particle potential is not 
exactly known, and polydispcrsity is inevitable. What 
is clear from the worked example summarised in Fig. [3] 
is that even in the 'ideal' case of monodispcrse spheres 
with an exactly-known inter-particle interaction, map- 
ping to hard spheres is system- and approach-specific. In 
practice, of course, polydispcrsity introduces significant 
uncertainties, and the Debye length is often not deter- 
minable to high accuracy. Moreover, we stress that the 
values of Z = 500 and Z ~ 100 are taken from experi- 
ments on nominally identical systems. Thus, conclusions 
derived from any 'mapping to hard spheres', e.g. com- 
parison of nucleation rates at nominally equivalent state 
points in the coexistence gap, must be treated with sig- 
nificant caution. 



VI. CONCLUSIONS 

We set out on a quest for colloids that mimic as closely 
as possible the ideal hard sphere. It seems that, to date, 
small (say a < 200nm) sterically PMMA particles dis- 
persed in index-matched hydrocarbons come about as 
close as possible to this ideal. Larger PMMA particles 
suitable for confocal microscopy require density match- 
ing to minimise sedimentation; but the solvents used 
inevitably induce a degree of charging that is difficult 
to screen out entirely using salts. The resulting soft, 
screened Coulomb inter-particlc repulsion can be satis- 
factorily modelled on the mean-field level, however the 
interaction parameters vary hugely for nominally identi- 
cal systems. The use of microgels such as PNiPAM, ei- 
ther on their own or as 'shells' on hard 'cores', has become 
popular, because their volume fraction can conveniently 
be tuned by temperature. Their inevitably soft mutual 
interaction has proven harder to model in a generic, an- 
alytical form. In both cases, the softness necessitates 
the use of a 4>oS to map onto hard spheres. We have 



reviewed various ways of performing this mapping, with 
or without the benefit of knowledge of the inter-particle 
potential u(r). None seem entirely satisfactory. 

We end by making two further observations. First, we 
widen the scope of our enquiry from charged particles 
and microgels to other kinds of non-hard inter-particle 
interaction, and ask what requirements should be satis- 
fied before one may fruitfully embark on the exercise of 
'mapping' to hard spheres. We suggest the the minimum 
conditions to be satisfied are: 

(1) the absence of any attractive interaction, so that the 
equilibrium physics is dominated by entropic effects; 

(2) crystals in a sufficiently monodisperse dispersion at 
high concentrations consist of the stacking of hexagonal 
layers. 

The first criterion highlights the importance of refractive 
index matching to minimise the ubiquitous vdW inter- 
action. The second criterion explains why c harg ed hard 
particles with no > 6, Fig. HI and microgels (54 1551] are 
suitable candidates for mapping to hard spheres, but star 
polymers are not.(70l| 

Secondly and finally, we point out that new develop- 
ments in particle synthesis may yet produce /xm-sized 
colloids that can be index and density matched using 
solvents or solvent mixtures that do not bring about 
charging and minimal swelling. Such development will 
be most welcome for the community of scientists wishing 
to use colloids to test fundamental theories of many-body 
physics via the hard-sphere model system. 

Of course, soft particles, from microgels to star poly- 
mers and beyond, are fascinating systems in their own 
right. [83| Furthermore, it is reasonable to enquire how 
much deviation from perfect hard spheres is acceptable, 
and the answer of course depends on what one wishes 
to study. Section [V] shows that the location of colloidal 
phase boundaries can vary in nontrivial and qualitative 
ways from hard spheres. However, it is plausible that 
slight softening may only have slight changes in, for ex- 
ample, the structural relaxation time [30| . 

Our point, then, is that while particle interaction de- 
tails may not matter in some cases (such as the pair struc- 
ture of dense liquids), there are plenty of cases where 
the behaviour of the system depends strongly on both 
the volume fraction cf> and the interparticle interactions, 
therefore accurate knowledge of the both is essential. 
Crucially, the interactions may well be known to even 
less precision that the absolute volume fraction <p, which 
we have already argued is knowable only to 3-6%. 
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